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ABSTRACT 

The Vela pulsar is the brightest pulsar at radio wavelengths. It was the 
object that told us (via its glitching) that pulsars were solid rotating bodies not 
oscillating ones. Along with the Crab pulsar is it the source of many of the 
models of pulsar behavior. Therefore it is of vital importance to know how far 
away it is, and its origin. 

The proper motion and parallax for the Vela pulsar have been derived from 
2.3 and 8.4 GHz Very Long Baseline Interferometry (VLBI) observations. The 
data spans 6.8 years and consists of eleven epochs. We find a proper motion of 
AW>s<5 = —49.61 ± 0.06, fis = 29.8 ±0.1 mas yr _1 and a parallax of 3.4 ± 0.2 
mas, which is equivalent to a distance of 293^ pc. When we subtract out the 
galactic rotation and solar peculiar velocity we find /z* = 45 ± 1.3 mas yr _1 with 
a position angle (PA) of 301° ±1.8 which implies that the proper motion has a 
small but significant offset from the X-ray nebula's symmetry axis. 



1. Introduction 

The Vela pulsar is an archetypal young pulsar. It was the first to be observed to glitch 
(Radhakrishnan & Manchester 1969) and was always associated with the large Vela su- 
pernova remnant (SNR). Since Vela is a young pulsar and displays a significant amount of 
glitching behavior and timing noise (Cordes et al. 1988), no accurate proper motion measure- 
ments from timing observations have been possible; it is only through direct proper motion 
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observations that we can trace its path on the sky. This is true of nearly all young pulsars, 
which are usually the most interesting. These have, for example, the possibility of discov- 
ering the associated birth SNR. The first accurate measurements of Vela's proper motion 
were produced using the Parkes-Tidbinbilla Interferometer (PTI) (Bailes et al. 1989). This 
275 km baseline at 1.6 GHz gave a resolution of 140 mas. This provided the first accurate 
measurement of the proper motion and thus, by extrapolation, the birth site. We have ex- 
tended this work to allow the third dimension to be derived via the parallax. Measurement 
of relative motion on milli-arcsecond scales is a challenging task and several attempts have 
been made, with various success. Previous observations have been made with ground based 
optical observations, with connected radio interferometers (PTI and the VLA) and the Hub- 
ble Space Telescope (HST). We now present results from the Australian Long baseline array 
(LBA). 

We have 18 observations of the position and can provide the most accurate measurements 
of parallax and proper motion. Recent measurements of the position of Vela against a 
number of background sources has been reported by De Luca et al. (2000); Caraveo et al. 
(2001). Our more accurate values, using completely independent methods, are compared 
with theirs, confirming the parallax values. Differences in the derivation of the proper motion 
are discussed. We find good agreement if we account for the probable galactic rotation of 
the reference stars. 

The distance to the Vela pulsar was originally estimated from the SNR distance by 
comparison to the Cygnus loop and IC443. This lead to an estimate of 500 pc (Milne 
1968), however it has long been argued that this is an overestimate. Our results have been 
foreshadowed by the predictions in several papers that have made the case for the distance to 
the Vela pulsar to be drastically reduced. Analysis of X-ray observations of the pulsar made 
with ROSAT (Page et al. 1996) and Chandra (Pavlov et al. 2001b), of the SNR also with 
ROSAT (Bocchino et al. 1999) and with optical absorption (Jenkins & Wallerstein 1995; Cha 
et al. 1999) have all suggested that the original distance estimate should perhaps be halved. 
All of these results, however, are to some extent model dependent and therefore doubt has 
existed over their accuracy. 

Furthermore, our results allow the testing of the prediction of Spruit & Phinney (1998) 
on the alignment of the spin axes of pulsars with their proper motion vector. The spin axis 
can be derived from the high resolution images from Chandra. These allow the symmetry, 
and thus presumably the spin, axes to be directly discerned. Only two cases have been 
tested so far, the Crab (Caraveo & Mignani 1999) and Vela. The extra accuracy we can 
provide refines the conclusion reached by Pavlov et al. (2000); Helfand et al. (2001) for the 
Vela alignment. A note of caution should be raised in that an alternate explanation for the 
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X-ray nebula around Vela has been put forward (Radhakrishnan & Deshpande 2001) which 
explains the structure in terms of the particle beam from the polar cap. This description 
requires that only the projected spin axis lies along the proper motion axis. If the model of 
Spruit & Phinney (1998) is correct any misalignment of the axes allows the us to estimate the 
timescale of the impulse in terms of the initial rotation period, a value of great importance 
to understanding the core collapse. The accuracy, therefore, of the alignment of the axes is 
of great theoretical importance. 



2. Observations and Data reduction 

The LBA array is an Australian national facility and is usually made up of six telescopes, 
three operated by the CSIRO (Australian Telescope Compact Array, Parkes and Mopra), two 
by the University of Tasmania (Hobart and Ceduna) and one in South Africa operated by the 
Hartebeesthoek Radio Astronomical observatory. In addition DSN telescopes at Tidbinbilla 
are often included. The observations reported here are between the Tidbinbilla 70m and 
the Hobart 26m antennae, a baseline of 832 km and at two 16-MHz wide frequencies; 2.29 
and 8.425-GHz. In addition there was one observation with the Australian array, 32-MHz 
centered on 8.417-GHz, although only the Hobart to Tidbinbilla data was directly used for 
the results in this paper. Details of all the telescopes and their parameters are available 
via the Australian Telescope National Facility (ATNF) website 4 . All data were recorded 
using beam- switching on a ~5 minute cycle, using the extra galactic phase-reference source 
Vela-G (a J20 oo = 08 h 33 m 22 s . 31563, 5 J20 oo = -44°41'38". 71463). Vela-G is part of the Radio 
Reference Frame and was an International Celestial Reference Frame (ICRF) candidate 
source (Ma et al. 1998). 

The first three epochs observations of the Vela pulsar (Table 1) were made using 
Mklll/MklllA VLBI recording systems, using recording mode B, which provides 14 con- 
tiguous frequency channels each 2 MHz wide in right-hand circular polarization at 2.3 GHz. 
Sampling was one-bit (two-level) for all data. None of the Mark III/IIIA observations used 
the pulsar gating nor binning described later. Only the Mklll positions were available for 
this analysis. The Mklll/MklllA VLBI data were processed at the Washington Mark IIIA 
correlator, located at the U.S. Naval Observatory (USNO). Phase models were applied with 
the CALC 6.0 package, supplying the best possible values of station locations, clock models 
and earth orientation parameters. The raw data were then exported in "FRNGX" format 
and further analyzed following the phase-referencing technique described by Lestrade et al. 
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(1990), and used the software SPRINT, developed by J.-F. Lestrade. 

The S2 recordings system and correlator for LB A were commissioned in the mid-1990s 
and observations were transferred to this system. In all observations, two separate frequency 
channels of 16 MHz bandwidth were recorded with two-bit (four-level) sampling. Most of 
the observations of the Vela pulsar using the S2 system were made with one 16 MHz channel 
at 2.3 GHz (13 cm) and one at 8.4 GHz (3 cm), taking advantage of the dual 13/3 cm 
(right-circular polarization) receivers at both Hobart and Tidbinbilla. In the final epoch, 
using the whole LBA at 8.4 GHz, two 16 MHz channels were used to form one contiguous 
32 MHz band. 

The S2 correlator is operated by ATNF at its Marsfield headquarters in Sydney. This 
wonderfully flexible device can pulsar bin (as opposed to pulsar gate) with thirty-two bins 
across the Vela pulsar's period, thereby limiting the dispersion to individual channels rather 
than the entire band. There is, however, no significant dispersion at these observing fre- 
quencies for Vela. We used only those pulsar bins containing significant flux, usually one or 
two out of thirty two, gaining up to a factor of 5.6 in signal to noise. 

Post correlation the data were fringe fitted at the nominal pulsar position using the 
binned flux, with ten minute solution interval in AIPS 5 , and then exported. Further pro- 
cessing was done with DIFMAP (Shepherd 1997). The data were flagged and averaged to 
generate a statistical weight. The offsets of Vela-G from its nominal reference position were 
found by fitting a point source to the visibilities. The position of Vela was shifted (phase 
rotated) to correct for these offsets. We then selfcalibrated the Vela phases in AIPS at 5 
minutes intervals, and copied those corrections to the Vela-G visibilities. The final offsets 
of Vela-G from the reference position, all very small, were used to calculate the final Vela 
position. 

Lobe ambiguity was not an issue even though this was a single baseline experiment, as 
observing with two widely separated frequencies broke any degeneracy. The final observation, 
which was with two 8.4-GHz bandpasses, had multiple baselines with which to identify the 
correct lobe. Only the Tidbinbilla to Hobart baseline was used for measurements of the 
positions. 
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3. Results 

We find a pulsar position of a j 2 ooo = 08 ft 35 m 20 s .61153±0.00002, <5 J20 oo = -45°10'34".8755± 
0.0003 for a reference epoch of 2000.0. The proper motion is /i aC os<5 = —49.61 ± 0.06, us = 
29.8 ± 0.1 mas yr _1 and a parallax of 3.4 ± 0.2 mas, equivalent to a distance of 293^ pc. 

3.1. Sources of error in phase referenced observations 

3.1.1. Source structure of the phase reference source 

Whilst strictly the pulsar itself is the phase-reference in this experiment (with a binned 
flux of 5 and 0.6 Jy at 2.3 and 8.4 GHz) the effect of significant structure is the same in 
either case. If the source is non-symmetrical the incoming wavefront arriving at the baseline 
will not have the assumed intrinsic zero phase, degrading the source fitting. Therefore for 
our final observation we used the full array to allow imaging of both Vela-G and the pulsar. 
As we were looking for local structure we were able to use selfcalibration. 

Following the Chandra observations which unveiled the fine X-ray structure around 
the pulsar we carefully searched for radio structure at VLBI resolutions in the off pulse 
data (calibrated with the on pulse phases). We found no discernible structure on baselines 
between 2.6 and 15.6 MA (equivalent to a full width half maximum of 16 - 95 mas) within 
10" of the pulsar. The peak point source flux was 0.8 mJy/beam in line with expectations 
given the RMS of 0.1 mJy This is of no great surprise as the radio structure around the 
pulsar has been reported to resolve out at resolutions finer than a few arcseconds (Dodson 
et al. 2003 Accepted; Lewis et al. 2002). 

As pulsar radiospheres have angular sizes much less than 1 mas, they may be treated 
as a point source. (Scattering caused by the ISM may cause observable angular broadening 
of the pulsar image but is unlikely to bias the observed position significantly). The phase- 
reference source has not been optically identified but is most likely a distant AGN, and might 
well have source structure that produces systematic errors in the positions determined at 
each epoch for the pulsar. The fact that in most observations we used the Vela pulsar as 
the reference source for the weaker Vela G does makes no essential difference in this respect. 
With this in mind, the final epoch of our observations was used to image the reference source 
with the full LBA array. No structure brighter than 1 mJy was found around Vela-G (18 
mJy) at 8.4 GHz, with an image RMS of 0.3 mJy. The phase residuals had a RMS of 4°. 
We note here that Vela-G is a well known extra-galactic source and considered for the ICRF 
(candidate source ICRF J083322.3-444138) (Ma et al. 1998). 
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3.1.2. Source structure of the ionospheric delay 

A considerable phase error can be caused by the different ionospheric delay encountered 
over the different antennae. Several approaches exist to address these problems, the best 
being using multiple frequencies in the VLBI observations and solving for the ionosphere 
as part of the reduction (Brisken et al. 2000). Unfortunately the S2 system used in our 
experiments does not have enough spanned bandwidth to allow this approach to be used. 
We therefore we have looked into using the measured Total Electron Count (TEC) from 
GPS observations (Walker & Chatterjee 1999). However the best data is on a 5° grid. The 
separation of our phase-reference and the target source is 0.7° and therefore GPS data cannot 
provide a useful correction. 

As no correction was possible we modeled the data quality using DIFWRAP (Lovell 
2000) which allows an error estimate that includes all possible contributions without at- 
tempting to identify them. This approach involves exploring a range of model parameters to 
identifying the range of acceptable fits, which should be the bounds of the multi-dimensional 
one sigma contour. The errors found are indeed approximately equivalent to one sigma, as 
confirmed by fact that the proper motion fit has a reduced \ 2 °f 11- 

3.2. Space velocity 

Two corrections need to be applied to our results to the true local motion of the Vela 
pulsar in its local environment. Our observations are directly tied to the ICRF, therefore the 
solar peculiar motion and the galactic rotation contribute to the observed proper motion, 
and need to be removed. We have used the solar constants from Dehnen & Binney (1998) 
of 10 ± 0.36,5.25 ± 0.62,7.17 ± 0.38 km s" 1 in galactic coordinates. We have used a flat 
rotation curve (Q = 220 km s -1 , R = 8.5 kpc; Fich et al. (1989)), which produces a local 
proper motion of —5.4 //; mas yr _1 . This compares with —5.7 \i\ mas yr _1 from the local 
values for the Oort constants, as found by Feast & Whitelock (1997). The corrections are 
shown in the Table 2. The dominant source of error is from the uncertainties in the solar 
peculiar motion parameters. We have used the measured uncertainties in our observations 
and combined those with the models. We have ignored the possibilities of systematic errors 
or alternate models. This gives us an angular motion, at Vela's local standard of rest, of 
/j, a = —38.5 ± 1.2 mas yr _1 , lis — 23 ± 1.5 mas yr^ 1 or //* = 45 ± 1.3 mas yr _1 at a position 
angle of 301° ± 1.8. The pulsars' transverse space velocity is therefore 62 ± 2 km s -1 . The 
PA no longer lies quite along the spin axis (e.g. 310° ± 1.5 Helfand et al. (2001), or 307° ±2 
Pavlov et al. (2001a)), which may strengthen the case that the impulse timescale was not 
quite long enough to average the off-axis component to zero. 
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4. Discussion 

4.1. Fitting methods 

We used the publicly available proper motion fitting routines, PMPAR 6 , created by W. 
Brisken. As each epoch had quite different observation spans and phase stability great care 
was required to ensure that the error estimates were accurate. The traditional approach has 
been to assume that the errors are a fraction of the beam size, but these ignore variations 
in the observing conditions during the experiment. We used DIFWRAP (Lovell 2000) to 
measure the complete range of errors and derived error estimates which are realistic. Where 
we only had the archival positions, and not the data, we have used the median value of the 
difference between the formal errors and the errors found with DIFWRAP. It is particularly 
important to get the errors correct, as the effect to be measured is small and the variation 
between the data quality is large. 

De Luca et al. (2000) point out that, as they calculate the proper motion using data 
collected on nearly identical day-numbers, their proper motion measurement is not contam- 
inated by the parallax. They refer to this as a 'pure' proper motion. This is a concern with 
very short time baselines, which would blend the parallax and proper motion. We, however, 
have data spanning seven years and the correlation between the parameters is low. 

4.2. Comparison of proper motion with other studies 

Our calculated proper motion and parallax can also be compared with recent values 
obtained by optical proper motion studies and phase-referencing VLBI. We ignore the historic 
observations, which were blighted by poor resolution and low elevation, and concentrate on 
the two Radio VLBI observations, this one and Bailes et al. (1989) and the four optical 
observations, one purely ground based (Nasuti et al. 1997), one ground and HST (Markwardt 
& Ogelman 1994) and two reports from the HST data set (De Luca et al. 2000; Caraveo 
et al. 2001) for which we take only the latest results. The results found by each of these 
studies is shown in Table 3. The table shows the proper motions in right ascension and 
declination for the pulsar along with the reference for this work. The proper motions listed 
do not account for the rotation of the Galaxy nor the peculiar motion of the Sun in the local 
galactic potential being only with reference to the calibrators used. The radio reference is 
an extra-galactic source, and the optical references are a significant number of field stars. 



6 http: / / nacho.princeton.edu/walterfb/pmpar/pmpar.html 



- 8- 



Both styles of observations are internally consistent to within two standard deviations (see 
Figure 2), with improvement of errors over time. With the reduction of errors, however, the 
radio and optical results are steadily becoming less compatible. 

The major advantage our observation have over the HST ones (other than the formal 
resolution being approximately one hundred times better) is that the reference is tied to 
a well defined reference frame. For the optical observations only field stars were available. 
The distance to the reference stars in the HST must be significantly greater than that of 
the pulsar, otherwise they would have had observable parallaxes themselves. Not knowing 
a distance we have assumed that they lie between two and ten kpc and calculated what the 
apparent proper motion would be for these limits using the standard flat rotation curve. We 
find that the galactic rotation contributes between -5 - -3 mas yr _1 to the proper motion 
in I. As all the sources would have proper motions within this range the scatter is within 
the HST errors of 1 mas. This contribution was not included in the calculations of Caraveo 
et al. (2001), and when it is the two sets of observations are consistent. Figure 2 includes 
a line representing the range that having the reference source at two to ten kpc would have 
contributed to our result. This means that the radio results are now consistent with all the 
optical results. This contribution, of course, needs to be removed to produce the correct 
space velocity and position angle for the pulsar. 

5. Conclusions 

We have measured the proper motion and parallax of the Vela pulsar to an unprece- 
dented accuracy (/J, acos s = —49.61 ±0.06, //$ = 29.8 ±0.1 mas yr _1 , it = 3.4 ±0.2 mas), and 
have been able to convert these back to the space velocity and position angle of the pulsar 
in its local environment with greater precision that previously possible (62 ± 2 km s _1 at 
301° ± 1.8), because of the unambiguity in the radio reference frame. This allows the pre- 
cise comparison of the Vela X-ray nebula symmetry axis and the proper motion of the Vela 
pulsar, opening insights into the timescale of the core collapse processes. 
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Fig. 1. — The residual offsets in the position for the Vela pulsar, after subtraction of the 
proper motion, in RA and Declination. Solutions for 13cm (triangles) and 3cm (boxes) 
observations are in shown with error bars in red (13cm) and green (3cm). The errors are 
the range of acceptable fits of the to model to the data, as described in the text, and are 
significantly smaller for the 3cm observations, as would be expected. 
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Fig. 2. — Vela proper motion determinations in RA and Declination (mas per year) with 
2-cr error ellipses. These are labelled with the reference name and row number from Table 
2. The effect of adding the correction for the galactic rotation of calibrators for distances 
between 2 and lOkpc is shown for the our result. 
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1999 Oct 26 
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8425 


1999 Oct 26 
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8417 


2000 Feb 21 
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Table 1: Observational details for the data presented here. The image residual after the 
subtraction of the point source is also quoted. fMk III observations. 



In galactic coordinates 
Observed proper motion —53.6 ± 0.08, —21.9 ± 0.06 

+Correction for the solar motion 6.7 ± 1.6, 4.9 ± 0.3 
+Correction for the galactic rotation 5.35 ± 0.3, 0.0 
Converted back to RA and Dec -38.5 ± 1.2, +23.1 ± 1.5 

Table 2: Corrections to the observed Vela proper motion. 



Table 3: Vela proper motion determinations. Historic and inaccurate observations have not 
been included. These are the proper motions as seen against the calibrators, i.e. the optical 
and radio results are aligned to different reference frames. 



No. 


Method 


pi a (mas yr -1 ) 


[i 6 (mas yr -1 ) 


Reference 


1 


Phase-referencing VLBI 


-48 ±4 


34 ±2 


Bailes et al. (1989) 


2 


Ground-based optical & HST 


-41 ±3 


26 ±3 


Markwardt & Ogclman (1994) 


3 


Ground-based optical 


-47 ±3 


22 ±3 


Nasuti ct al. (1997) 


4 


HST 


-45 ± 1 


26 ± 1 


Caravco ct al. (2001) 


5 


This Result 


-49.61 ±0.06 


29.8 ± 0.1 





